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Foreword 

The  Wood-Witt  program  was  organized  in  1999,  and  has  involved  more  than  fifty  participants,  firom  more 
than  ten  different  countries.  The  present  grant,  F49620-03-1-0197,  has  covered  the  period  lMay03  - 
30Apr06.  The  overall  goal  has  been  to  understand  defects  and  impurities  in  GaN,  and  determine  how  to 
reduce  them,  if  possible.  More  than  600  GaN  samples  have  been  investigated  in  our  laboratory  or  sent  out 
since  the  inception  of  the  program,  and  much  has  been  learned  about  this  important  material.  The 
participants  from  Wright  State  University  (WSU)  during  the  last  three  years  have  included  Bruce  (Chip) 
Claflin,  Timothy  Cooper,  Zhaoqiang  Fang,  Gary  Farlow,  John  Hoelscher,  Brahmanand  (Val)  Jogai,  David 
Look,  Donald  Reynolds,  and  Wallace  Rice,  Jr.  There  have  also  been  participants  from  the  Materials  and 
Manufacturing  Directorate,  and  Sensors  Directorate  of  the  Air  Force  Research  Laboratory  (AFRL),  and 
many  others  from  around  the  world.  We  especially  want  to  thank  our  AFOSR  program  managers,  Dr. 
Gerald  L.  Witt,  who  initiated  the  program.  Dr.  Kitt  C.  Reinhardt,  who  took  over  near  the  end.  Their  insight 
and  support  has  been  crucial  in  making  this  program  successful.  One  measure  of  success  is  the  number  and 
quality  of  journal  publications  resulting  from  the  research.  In  the  last  three  years,  our  group  at  WSU  has 
generated  twenty-five  papers,  eo-authored  by  scientists  in  ten  different  countries.  Also,  many  additional 
papers  from  other  groups  have  been  published  on  the  same  set  of  samples.  A  list  of  the  WSU  publications  in 
included  at  the  end  of  this  report.  We  have  also  chosen  three  topics  to  illustrate  in  more  detail  in  the  body  of 
the  report,  below.  Each  of  the  three  sections  is  self  contained,  and  thus  can  be  read  without  reference  to 
either  of  the  other  two.  Any  questions  or  comments  on  the  research  can  be  directed  to  Dr.  David  Look, 
david.look@wright.edu. 
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The  nitrogen  vacancy  in  GaN 


1.1  Introduction 

The  last  decade  has  seen  greatly  increased  research  and  development  on  GaN-related  materials  and 
devices.'  The  driving  force  has  been  the  potential  for  blue/UV  light-emitting  diodes  (LEDs)  and  laser  diodes 
(LDs),  and  also  high-frequency  transistors  operating  at  high  powers  and  temperatures.  Much  of  the  GaN 
research  has  concentrated  on  identitying  impurities  and  defects  that  may  act  as  donors,  acceptors,  traps,  or 
recombination  centers.  In  particular,  the  main  donors  that  have  been  investigated  so  far  are  oxygen 
occupying  a  nitrogen  site  (On),  silicon  on  a  Ga  site  (Sica),  and  the  N  vacancy.^’''  Several  years  ago,  we 
performed  0.7  -1.0  MeV  electron  irradiation  on  the  best  (highest-mobility)  GaN  samples  available  at  that 
time,  and  found  that  a  donor  (at  64  ±  10  meV  below  the  conduction-band  edge)  and  an  acceptor  (much 
lower  in  the  bandgap)  were  produced  at  approximately  the  same  rate.  We  then  used  several  arguments  to 
assign  the  donor  to  the  N  vacancy,  and  the  acceptor  to  the  N  interstitial.''  Recently,  GaN  samples  with  much 
lower  background  donor  and  acceptor  concentrations  have  become  available,  allowing  better  accuracy  in  the 
Hall  fitting,  and  sharper  and  more  intense  photoluminescence  (PL)  spectra.^'^  Furthermore,  we  have  gained 
access  to  a  lower-energy  electron  accelerator,  allowing  a  separation  of  N-sublattice  from  Ga-sublattice 
damage.  Thus,  we  have  revisited  the  GaN  defect  problem  by  irradiating  this  new,  higher-quality  GaN  with 
0.42-MeV  electrons,  which  recent  displacement-energy  calculations  have  shown  to  be  above  the  N- 

o 

sublattice  damage  threshold,  but  below  that  of  the  Ga  sublattice.  Therefore,  in  this  case,  we  can  be  sure  that 
only  N  vacancies  (Vn)  and  interstitials  (Ni)  are  initially  being  produced  by  the  irradiation.  However,  the  Ni 
are  likely  mobile  at  room  temperature,  either  recombining  with  vacancies,  or  forming  complexes  with 
existing  impurities  or  defects.  One  of  these  irradiated  samples  was  given  to  a  group  performing  PL 
measurements,  and  new,  weak  lines  were  seen  in  the  donor-bound  exciton  (DBE)  region  and  the  two- 
electron  satellite  (TES)  region.  By  assuming  that  the  DBE  and  TES  lines  represented,  respectively,  n=l  and 
n=2  final  states  of  the  same  donor,  it  was  concluded  that  a  25-meV  donor  was  being  created  by  the 
irradiation.  Furthermore,  this  donor  was  assigned  to  the  N  vacancy.  However,  we  show  below  that  this 
assignment  is  completely  incompatible  with  our  Hall  results,  since  donors  in  the  25-meV  region  actually 
decrease  in  concentration,  and  is  also  incompatible  with  our  PL  results,  since  no  new  DBE  lines  are 
observed  in  our  irradiated  sample. 

1.2  Experimental  considerations 
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The  GaN  sample  diseussed  here  was  grown  in  the  (0001)  orientation  (Ga  face  up)  on  AI2O3  by  the 
hydride  vapor-phase  epitaxial  (HVPE)  technique  at  the  Samsung  Advanced  Institute  of  Technology.^ 
Separation  of  the  GaN  and  AI2O3  was  effected  by  laser  irradiation  on  the  N  face,  through  the  AI2O3 
substrate.  Lapping  and  polishing  of  both  the  Ga  and  N  faces  produced  a  flat  sample  of  thickness  219  pm. 
This  process  also  removed  the  highly-conductive  N-face/Al203  interface  layer,  which  made  analysis 
difficult  in  our  first  study Samples  of  this  type  have  been  extensively  characterized  by  optical,  electrical, 
and  structural  techniques,  and  have  demonstrated  record  mobilities. 

The  van  der  Pauw  Hall-effect  measurements  were  performed  with  a  LakeShore  Model  7507  apparatus, 
including  a  closed-cycle  He  cooling  system  operating  from  15  -  320  K.  From  measurements  of  Hall 
coefficient  R  and  conductivity  a,  the  Hall  mobility  ph  =  Re  and  the  Hall  concentration  nn  =  1/eR  could  be 
calculated  at  each  temperature.  The  true  carrier  concentration  n  is  related  to  nn  by  n  =  mn,  where  r  is  the  so- 
called  Hall  r-factor.^  For  the  unirradiated  sample,  r  =  1.1  -  1.3,  from  50  -  320  K,  and  for  the  irradiated 
sample,  r=  1.1  -  1.5. 

Photoluminescence  measurements  were  performed  at  4.2  K.  Excitation,  dispersion,  and  detection  were 
accomplished,  respectively,  with  a  45-mW  HeCd  laser,  a  1.25-m  spectrometer,  and  a  photomultiplier 
detector.  Resolution  was  better  than  0.01  meV  in  the  spectral  range  important  for  this  study.  Electron 
irradiations  were  carried  out  at  room  temperature  with  the  beam  directed  in  the  [000-1]  direction,  i.e., 
opposite  to  the  growth  direction.  The  energy  was  about  0.42  MeV,  the  current  about  2  pA/cm  ,  and  the  total 
fluence  about  3.6  x  lO'^  cm'^.  Very  recently,  threshold  displacement  energies  (Ed’s)  have  been  calculated 
for  GaN  using  realistic  potentials.*  Depending  on  irradiation  direction,  the  minimum  Ed  for  N  displacement 
is  25  eV,  and  that  for  Ga  displacement,  22  eV.  However,  along  the  [000-1]  direction,  and  averaging  over  a 
15  °  acceptance  angle  to  account  for  thermal  motions  and  possible  beam  misalignment,  the  calculated  Ed’s 
are  66  and  38  eV,  respectively,  for  N  and  Ga  displacements.''^  From  these  theoretical  results,  and 
accounting  for  the  higher  mass  of  Ga,  it  is  found  that  the  minimum  electron  energy  needed  for  N 
displacement  is  0.32  MeV,  and  for  Ga  displacement,  0.53  MeV.  Furthermore,  using  the  McKinley- 
Feshbach  relativistic  displacement-cross-section  formula,'*  the  N  production  rate  should  be  about  0.03  cm’' 
at  our  electron  energy  of  0.42  MeV,  and  will  exceed  the  Ga  production  rate  up  to  electron  energies  of  0.87 
MeV.  These  considerations  support  our  conclusion  that  the  present  irradiations  are  primarily  producing  N- 
sublattice  displacements.  Ga-sublattice  displacements  have  been  identified  in  irradiations  at  much  higher 
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energies,  e.g.,  2.0  -  2.5  MeV.  ’  However,  the  Ga  interstitial  Gai  is  mobile  at  room  temperature,  and  tends 
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to  form  complexes.’^  Furthermore,  only  tightly-bound  wave  functions  have  been  associated  with  Gai."’’^ 
Thus,  the  isolated  Gai  is  not  a  good  candidate  for  the  shallow  donor  produced  by  room-temperature 
irradiation,  and  the  Ga  vacancy  is,  of  course,  not  a  donor  at  all,  but  a  well-known  acceptor.  From  these 
results,  only  the  N  vacancy  is  a  reasonable  candidate  for  the  shallow  donor  produced  by  0.42-MeV  electron 
irradiation. 

1.3  Results 

The  temperature  dependence  of  carrier  concentration  n  and  mobility  p  are  shown  in  Fig.  1.  The  peak 
mobility  is  mainly  sensitive  to  the  acceptor  concentration  Na,  showing  that  Na  increases  with  irradiation.'^’^ 
Note  that  Na  is  the  total  acceptor  concentration,  including  the  original  acceptors  as  well  as  those  produced 
by  the  irradiation.  (From  the  Hall-effect  measurements,  we  cannot  determine  the  energy  levels  of  the 
acceptors,  except  that  they  must  be  at  least  several  kT  below  the  donor  levels.)  In  Fig.  1,  it  is  clearly  seen 
that  another,  deeper  donor  Nd2  has  been  produced  by  the  irradiation,  although  the  original,  shallower  donor 
Ndi  is  also  still  present.  The  fitted  values  of  Not,  Edi,  Ndi,  Ed2,  and  Na  are  given  in  Table  1.  Here,  Na  is 
determined  from  the  mobility  data,  and  all  the  other  parameters  from  the  carrier-concentration  data.^  Note 
that  a  donor  at  70  ±  2  meV  increases  in  concentration  by  7.0  x  lO’^  cm"^,  giving  a  production  rate  of  about 
0.02  cm"',  close  to  the  predicted  value  of  0.03  cm"'.  In  contrast,  the  donors  at  25  ±  1  meV  actually  decrease 
in  concentration,  by  1.8  x  lO'^  cm‘^.  This  observation  can  be  understood  by  realizing  that  the  N  interstitials 
are  likely  mobile  at  room  temperature,  and  will  tend  to  form  complexes  with  impurities  or  other  defects. 
Since  the  main  donor  impurity  in  the  present  material  is  probably  On,^  a  single  donor,  and  since  Ni  is 
thought  to  be  a  single  acceptor  in  n-type  GaN,’'*  there  would  be  a  coulomb  attraction  between  On  and  Ni, 
and  the  resulting  complex  On-Ni  would  probably  be  neutral.  This  process  could  explain  why  the  25-meV 
donors  decrease  in  concentration.  In  this  scenario,  the  remaining  Ni  (of  concentration  5.2  x  lO'^  cm'^) 
would  account  for  the  increase  in  acceptor  concentration  (3.0  x  lO'^  cm'^),  and  also  a  possible  increase  in 
some  other  neutral  centers,  not  associated  with  the  original  shallow  donors. 

It  is  interesting  to  compare  the  present  donor  energy  of  70  ±  2  meV  with  the  value  64  ±  10  meV 
dedueed  in  the  earlier  irradiation  study.'*  In  the  earlier  case,  the  background  donor  concentration  was  1.2  x 
lO'^  cm’^,  whereas  the  background  donor  concentration  in  the  present  sample  is  1.25  x  lO’^  cm'^.  Using  the 
screening  formula,  Ed  =  Edo  -  aND*^^,  with  a  =  2.1  x  10'^  meV-cm,'^  we  get  an  unsereened  energy  Edo  =  75 
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meV  in  the  present  study,  and  74  meV  in  the  former  study.  Thus,  there  is  no  doubt  that  the  defect  donors 
produced  at  0.42  MeV,  and  those  produced  at  0.7  -  1 .0  MeV,  are  the  same. 

1.4  Discussion 

A  recent  paper  by  Yang  et  al.  (henceforth  called  YFW),'^  based  on  a  nearly  identical  irradiated 
sample,  concludes  that  an  optically-observed  25-meV  center  is  the  N  vacancy.  This  assignment  is 
immediately  incompatible  with  our  results  in  that  the  donors  at  25  +  1  meV  actually  decrease  in 
concentration.  Furthermore,  we  can  compare  the  YFW  PL  data  with  our  own  PL  data.  In  Fig.  2,  we  present 
donor-bound  exciton  spectra  from  three  different  samples:  (1)  unirradiated  S422  (dashed  line);  (2)  irradiated 
S422  (solid  line);  and  (3)  unirradiated  S417  (chained  line).  Here,  sample  S417  is  shown  for  comparison, 
since  YFW  had  an  adjacent  piece  of  S417.  Note  the  three  well-defined  DBE  lines  in  S417,  at  3.47136, 
3.47241,  and  3.47319  eV,  respectively.  The  identity  of  the  first  of  these  is  unknown,  whereas  the  other  two 
are  usually  assigned  to  On  and  Sioa,  respectively.^’’*'’  The  unirradiated  and  irradiated  S422  samples  exhibit 
the  same  three  lines,  but  they  are  not  as  sharply  defined.  In  their  unirradiated  sample  S417,  YFW  find  the 
dominant  line  (On)  at  3.4717  eV,  and  the  Sioa  line  at  3.4725  eV,  each  about  0.7  meV  below  our 
corresponding  values,  due  to  a  different  spectrometer  calibration.  But  after  irradiation,  they  find  a  new  line, 
at  3.4732  eV.  The  On,  Sioa,  and  new  line,  each  raised  by  0.7  meV  to  match  our  On  and  Sioa  lines,  are 
schematically  represented  by  solid  vertical  lines  in  Fig.  2.  Clearly,  we  do  not  see  their  new  DBE  line  in  our 
irradiated  sample.  However,  they  also  observe  another  weak  line  produced  by  the  irradiation,  at  3.4547  eV 
(not  shown),  and  this  they  believe  is  the  so-called  two-electron  satellite  (TES)  line  of  the  new  DBE  line  at 
3.4732.  (Note  that  TES  lines  result  from  bound-exciton  transitions  in  which  the  donor  is  left  in  an  n=2  state.) 
Then,  by  applying  the  hydrogenic  model  to  the  DBE  (n=l)  and  TES  (n=2)  lines,  the  defect  donor  energy 
becomes  4/3(3.4732  -  3.4547)  =  25  meV.  Indeed,  we  also  see  a  very  weak  feature  at  about  3.456  eV  (not 
shown),  perhaps  corresponding  to  their  feature  at  3.4547  eV.  However,  because  we  see  no  corresponding 
DBE  (n=l)  line,  we  cannot  interpret  this  feature  as  a  TES  (n=2)  line.  Instead,  it  perhaps  represents  a  donor- 
electron  to  free-hole  transition,  with  the  donor  having  an  energy  of  about  50  meV.  If  it  indeed  exists,  this  50- 
meV  donor  could  arise  from  complexes  formed  from  the  migrating  Ni  defects.  However,  we  reiterate  that 
neither  this  possible  center,  nor  the  25-meV  center  postulated  by  YFW,  have  concentrations  high  enough  to 
be  observed  by  the  Hall-effect  measurements.  In  contrast,  the  70-meV  donor  is  strongly  produced,  with  a 
concentration  that  is  close  to  the  theoretically  predicted  one;  thus,  only  the  70-meV  donor  is  a  good 
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candidate  for  the  isolated  N  vacancy.  In  a  further  comment  on  PL  data  in  irradiated  GaN,  we  may  note  that 
if  the  dominant  On  DBE  line  is  reduced  by  pairing  with  migrating  Ni  atoms,  then  the  nearby  Sica  DBE  line 
and  possibly  other  weaker  DBE  lines  may  appear  to  be  relatively  larger,  as  seen  in  our  PL  data  (Fig.  2). 
Depending  on  spectrometer  resolution  and  other  factors,  it  may  be  mistakenly  assumed  that  new  lines  are 
appearing. 

A  final  argument  by  YFW  is  that  their  25-meV  center  is  a  good  candidate  for  the  N  vacaney,  because 
this  vacancy  would  be  expected  to  have  an  energy  less  than  that  of  On,  whieh  is  known  to  be  nearly 
hydrogenic  with  an  energy  of  about  33  -  34  meV.  ’  (Note  that  the  ~25-meV  donor  measured  by  the  Hall 
effect  is  equivalent  to  the  ~33-meV  donor  measured  by  PL,  due  to  wavefunction  overlap  and  other  factors 
which  effectively  reduce  the  Hall  energy.^  )  That  is,  they  have  suggested  that  the  positive  charge  of  the  N 
vacancy  would  be  distributed  over  the  electron  states  of  the  surrounding  Ga  atoms,  and  thus  the  donor 
electron  would  require  less  energy  for  ionization.  However,  this  simple  argument  is  strongly  violated  in  a 
very  similar  case,  that  of  the  As  vacancy  in  GaAs.  In  GaAs,  the  hydrogenie  (0/+)  donor  energy  is  about  6 
meV,  whereas  the  (0/+)  transition  of  the  As  vacaney  lies  at  about  140  meV.’*’^°  Even  the  (-/O)  acceptor 
transition  of  the  As  vacancy,  at  45  meV,  is  deeper  than  the  hydrogenic  donor  level  in  GaAs.  From  this 
example,  the  N  vacancy  in  GaN  would  be  expected  to  be  deeper  than  the  hydrogenie  donor  level  (as  we 
observe),  not  shallower.  However,  further  consideration  of  the  Vn  energy  must  await  accurate  theoretical 
modeling. 
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Table  I.  Fitted  concentrations  and  energies  for  sample  S422,  irradiated  and  unirradiated. 


Irradiation 

Ndi  (cm’^) 

Edi  (meV) 

ND2(cm-®) 

ED2(meV) 

NAlcm"®) 

none 

1.25x10^® 

25.9 

1.19x10^® 

72.2 

2.30x10^® 

0.42  MeV 

1.07x10^® 

23.9 

8.18x10^® 

68.5 

5.29x10^® 
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Fig.  1.  Carrier  concentration,  corrected  for  Hall  r-factor,  for  unirradiated  and  irradiated  sample  S422.  Inset: 
Hall  mobility,  in  units  of  cm^A^-s.  All  solid  lines  are  theoretical  fits  to  the  experimental  data  (points). 
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Fig.  2.  Photoluminescence  donor-bound  exciton  region  for  sample  S422  (unirradiated  and  irradiated) 
and  sample  S417  (unirradiated).  Intensities  are  normalized  to  the  peak  intensity  of  the  3.47241-eV  line. 
The  three,  short  vertical  segments  at  the  bottom  of  the  figure  represent  lines  observed  by  the  authors  of 
Ref  15  (YFW).  The  two  lowest-energy  segments  represent  PL  lines  in  their  unirradiated  sample, 
whereas  the  highest-energy  segment  represents  a  line  observed  by  them  only  after  irradiation. 
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2.0  Giant  traps  associated  with  extended  defects  in  GaN  and  SiC 
2.1  Introduction 

It  is  well-known  that  impurities  and  point  defeets  can  act  as  donors,  acceptors,  recombination 
centers,  and  traps  in  semiconductor  materials.  For  simple  point  charges,  models  to  describe  their  effects  on 
electrical  and  optical  properties  are  readily  available.  Traps,  the  subject  of  this  paper,  are  often  investigated 
by  deep  level  transient  spectroscopy  (DLTS),  and  the  usual  analysis  of  DLTS  data  is  based  on  the 
assumption  that  each  trap  has  a  constant  capture  cross  section,  independent  of  time  or  the  presence  of  nearby 
traps  [1-3].  Thus,  if  the  Fermi  level  is  suddenly  raised,  say  by  applying  a  forward  bias  to  a  Schottky  contact, 
electron  traps  will  capture  free  electrons  in  an  exponential  manner,  and  if  the  Fermi  level  is  then  returned  to 
its  original  value,  the  traps  will  emit  the  electrons  in  the  same  way.  However,  large  defects,  such  as 
dislocations,  can  often  hold  many  charges,  and  they  are  not  independent  of  each  other  [4-8].  In  the  case  of  a 
threading  dislocation,  the  charges  can  build  up  along  a  line,  and  gradually  begin  to  repel  the  accumulation  of 
additional  charges.  Thus,  the  effective  capture  cross  section  decreases  as  the  accumulated  fraction  becomes 
larger,  and  the  capture  process  is  no  longer  exponential  with  time.  The  dynamics  of  dislocation  capture  were 
first  elucidated  in  1989  [5],  and  since  then  have  been  applied  to  many  experimental  situations.  Recently,  it 
has  been  shown  that  pores  also  can  capture  charge  [7,8],  and  that  both  the  capture  and  emission  processes 
are  quite  nonexponential  [8].  To  model  this  case,  a  new  formalism  has  been  developed  [8],  and  it  is 
applicable  to  any  charge  configuration  that  can  be  approximated  by  a  sphere,  cylinder,  or  thin  slab. 
Examples  of  such  shapes  might  include  pores,  open-core  dislocations,  and  cracks,  respectively.  In  this  work, 
we  specifically  investigate  pores  in  SiC,  and  dislocations  and  microcracks  in  GaN. 

2.2  Theory 

Consider  a  spherical  pore  of  radius  rp  with  deep,  single-level  acceptor  states  of  sheet  density  Nss  on 
its  inner  surface.  If  all  of  the  states  are  filled,  then  the  total  number  of  traps  is  4Tcrp^Nss.  However,  as  more 
and  more  electrons  are  trapped,  a  negative  (repulsive)  potential  (jisph  builds  up,  and  the  trapping  rate 
diminishes.  A  spherical  region  depleted  of  free  electrons,  described  by  a  local  band  bending  of  energy  Oisph 
6<t^sph(rp),  forms  outside  the  surface  (r  =  rp)  of  the  pore.  The  value  of  O  can  be  calculated  from  Poisson’s 
equation,  which,  for  spherical  pores,  is  most  conveniently  expressed  in  spherical  coordinates: 
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where  p  is  the  charge  density,  e  is  the  static  dielectric  constant,  and  No  is  the  net  donor  density  (actually.  No 
-  Na,  where  Na  is  the  acceptor  concentration).  By  applying  charge  conservation  and  the  depletion 
approximation,  it  can  be  shown  that  the  energy,  O  =  -e(|),  is  given  by 
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where  f  is  the  fractional  occupation  of  the  trap  states  on  the  pore,  i.e.,  f  =  Nss’/Nss-  Equation  2  holds  for  rp  <  r 
<  w,  where  w  is  the  depletion  length,  and  is  cast  in  a  form  which  is  convenient  in  that  the  first  two  terms 
drop  out  for  r  =  tp.  (For  r  >  w,  O  =  0.)  For  our  porous  SiC  sample.  No  »  lO'^  cm'^  and  rp «  20  nm,  and  it  will 
turn  out  that  Nss »  2.5  x  lO'^  cm'^,  so  that  0(rp) «  0.2  eV,  for  f  =  1 . 

For  cylindrical  pores,  it  can  be  shown  that 
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for  rp  <  r  <  w.  In  this  case,  0(rp) «  0.3  eV,  for  f  =  1 . 

The  dynamic  capture  and  emission  processes  can  now  be  described  by  the  usual  master  equation: 

^  =  +  (4) 

where  Cn  is  the  emission  rate  from  filled  traps,  and  Cn  is  the  capture  rate  into  empty  traps.  Both  rates  will  be 
affected  (slowed  down)  by  the  band  bending,  as  follows: 
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where  v(T)  =  (8kT/7im*)''^^  is  the  thermal  velocity,  Ncb(T)  =  2(27Tm*kT)^^^/h^  is  the  effective  conduction 
band  density  of  states  (in  the  Boltzmann  approximation),  n  is  the  free  electron  concentration,  a  is  the 
capture  cross  section  for  a  single  trap,  and  Gemis  -  (go/gi)<7exp(a/k),  where  go  and  gi  are  the  degeneracies  of 
the  unoccupied  and  occupied  trap  states,  respectively,  and  a  is  a  linear  temperature  coefficient:  Ess  =  Esso(l 
-  aT).  From  Eqns.  4  -  6,  a  general  integral  equation,  describing  both  capture  and  emission,  can  be  written 
for  f(t): 
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In  a  DLTS  experiment,  the  capture  process  is  carried  out  by  applying  a  forward  bias  to  a  normally 
reversed-biased  Schottky  barrier  or  p-n  junction  [1-3].  In  reverse  bias,  the  traps  are  in  a  region  depleted  of 
free  electrons,  and  thus  experience  a  very  low  free-electron  concentration,  n  =  nr «  nb,  where  nb  is  the  bulk 
(neutral)  value,  lO'*  cm‘^  in  this  case.  Thus,  Cno(nr,T)  is  very  small,  so  that  emission  dominates  and  the  traps 
are  almost  empty.  Then,  in  forward  bias,  the  traps  are  suddenly  exposed  to  the  bulk  free-electron 
concentration  n  =  nb  for  a  time  tp,  the  filling  pulse  length,  and  at  the  end  of  this  pulse  the  filled  fraction  is 
defined  as  fp.  Thus,  the  trap  filling  process  is  described  by  solving  Eq.  7  for  fp  under  the  conditions  fa  =  0,  fp 
=  fp,  n  =  nb,  ta  =  0,  and  tp  =  tp.  (One  convenient  means  of  solving  Eq.  7  is  by  use  of  the  “root”  function  in 
Mathcad  [9].) 

When  the  filling  pulse  has  ended  at  time  tp,  i.e.,  by  reapplying  the  reverse  bias,  the  traps  are  once  again 
suddenly  exposed  to  a  very  small  value  of  n,  i.e.,  n  =  nr.  (Note  that  the  solution  of  Eq.  7  is  very  insensitive  to 
the  exact  value  of  nr,  as  long  as  nr  «  nb.)  The  traps  now  emit  their  carriers,  so  that  the  original  fractional 
occupation  fp  is  now  reduced  to  fe,  in  total  time  tp  +  U-  Thus,  in  emission,  Eq.  7  is  solved  for  fe  under  the 
conditions  fa  =  fp,  fp  =  fe,  n  =  nr,  ta  =  tp,  and  tp  =  tp  +  te.  In  the  most  common  form  of  DLTS  methodology. 
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used  in  commercial  instruments  and  often  called  the  “boxcar”  technique  [1,2],  the  emission  curve  is 
evaluated  at  two  points,  ti  and  t2  (referenced  to  tp),  and  the  signal  strength  is  measured  as  S  =  f(ti)  -  f(t2). 
Such  a  signal  is  simulated  simply  by  solving  Eq.  7  at  two  times,  tp  +  ti,  and  tp  +  t2. 

To  illustrate  these  ideas,  we  show  in  Fig.  1  a  band  profile  along  a  line  perpendicular  to  the  sample 
surface  and  cutting  through  a  single,  spherical  pore,  of  radius  2  nm,  and  of  depth  22  nm  (from  the  surface). 
Here  the  Fermi  energy  Ef  in  the  bulk  (  depth  z  =  oo)  is  set  at  zero,  and  in  this  region  Ef  is  very  close  to  the 
conduction-band  energy  Ecb,  at  least  for  this  particular  sample.  While  under  reverse  bias  Vr  =  -V,  the 
surface  will  be  held  at  energy  e((t)B  -  Vr),  where  (|)b  is  the  Schottky-barrier  potential,  and  Ecb  will  vary 
smoothly  from  e((j)B  -  Vr)  at  z  =  0,  to  Ecb  =  0  at  z  >  w,  the  depletion  depth.  During  the  filling  pulse,  it  is 
assumed  (in  this  case)  that  Ecb  =  0  all  the  way  to  the  surface  (flat-band  approximation),  although  this 
assumption  is  not  invoked  for  the  actual  calculations.  The  dashed  lines  in  the  figure  illustrate  the  potential 
barrier  formed  near  the  surface  of  a  pore  due  to  the  captured  charge.  This  barrier  begins  to  build  up  as  soon 
as  the  filling  pulse  is  applied,  and  then  recedes  when  the  filling  pulse  is  removed.  The  changing  barrier 
energy  means  that  the  effective  capture  and  emission  rates  will  change  continuously  with  time,  as  described 
by  Eq.  7. 

Before  going  on,  it  is  instructive  to  solve  Eq.  7  for  a  particularly  simple  case,  0(rp,f)  =  0  (or  a  constant), 
which  should  hold  for  isolated  point  charges.  The  Eq.  7  yields  the  usual  closed-form,  exponential  solutions 
for  capture  and  emission: 


fcap,iT,tJ  =  - 


'capl\  ’-p-f  g  ^2"^ 


l-e 


-K(r)+c, 
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(9) 


These  equations  form  the  basis  of  the  original  DLTS  analysis  method  [1,2],  and  are  still  used  in  the  vast 
majority  of  cases,  as  well  as  in  commereial  instruments.  Besides  isolated  point  charges,  line  charges  have 
also  been  considered  in  the  literature  [5,6],  and  the  derived  capture  process  in  this  case  turns  out  to  be 
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logarithmic  with  time.  However,  such  a  time  dependence  is  only  an  approximation,  and  to  derive  it  from  the 
general  formalism,  Eq.  7,  we  must  assume  (as  did  Wosinski^)  the  following:  (1)  small  f,  such  that  the 
denominator  of  the  integrand  in  Eq.  7  can  be  approximated  by  unity;  and  (2)  0(rp,f)  oc  f,  which  means  that 
the  large,  bracketed  terms  in  Eqns.  2  or  3,  respectively,  must  be  independent  of  f,  i.e.,  a  constant  K.  Then, 
Eq.  7  yields  a  logarithmic  solution  for  f  (cf  Eq.  3  of  Ref  6),  which  has  been  seen  experimentally  for 
trapping  along  dislocation  lines [5, 6]: 
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where  Ksph  is  typically  on  the  order  of  0. 1  -  1 . 

2.3  Pores  in  SiC 

The  porous  SiC  in  this  study  was  prepared  by  photo-assisted  electrochemical  etching  [10]  of  n-type  6H  SiC 
obtained  from  Sterling  Semiconductor  Inc.  (now  part  of  Dow-Coming  Corp.).  The  electrolyte  was  a  mixture 
of  HE  acid  and  ethanol.  The  resistivity  of  the  starting  material  was  about  0.2  Q-cm,  and  the  carrier 
concentration  ~  lO'*  cm'^.  The  C-V  and  DLTS  data  were  obtained  by  means  of  a  BioRad  DL4600  DLTS 
apparatus,  which  operated  over  the  temperature  range  80  -  450  K.  From  the  C-V  data,  the  carrier 
concentration  in  the  nonporous  SiC  was  uniform  at  about  10  cm' ,  whereas  that  in  the  porous  SiC  dropped 
to  about  10  cm'  at  a  depth  of  about  80  nm.  Cross-sectional  transmission  electron  microscopy  (TEM), 
using  a  200kV  Phillips  CM-200  insfrument,  was  used  to  study  the  pore  size  and  density.  At  depths  of  50  - 
100  nm  below  the  surface,  the  pore  radii  ranged  in  size  from  10-25  nm,  with  a  density  of  about  5  x  lO'^ 
cm'^,  increasing  with  depth.  The  sizes  and  densities  of  these  pores  are  very  typical  of  those  found  at  the 
same  depth  (just  below  the  so-called  “skin  layer”)  in  other  P-SiC  samples  [7,10,1 1]. 

In  Fig.  2,  we  compare  the  capture  and  emission  solutions  under  three  different  approximations:  (1)  exact 
analysis  (Eq.  7);  (2)  exponential  analysis  (setting  Osph  =  0,  in  Eq.  7);  and  (3)  logarithmic  analysis  (for  f «  1, 
and  Osph  =  Ke  Nssrpfrs  in  Eq.  7).  To  generate  the  curves,  we  have  used  some  SiC  parameters  from  the 
literature:  m*/mo  =  0.4,  and  e/sq  =  10;  some  parameters  measured  by  TEM  or  C-V:  rp  =  2  x  10'®  cm'^;  No  « 
nb «  10  cm' ;  and  rie «  10  cm'  (fit  not  sensitive  to  ng);  and  some  fitted  parameters  (i.e.,  those  needed  to  fit 
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the  DLTS  data  of  trap  To  in  Fig.  3):  Nss  =  2.5  x  cm'^;  Esso  -  0.8  eV;  a  =  1  x  10'^^  cm^;  and  Oemis  =  3  x 
10’’^  cm^.  A  filling  pulse  length  tp  =  20  ms  was  assumed  for  the  curves  in  Fig.  3.  The  exponential 
approximation  is  the  one  assumed  in  the  vast  majority  of  DLTS  experiments,  and  indeed,  it  works  well  for 
simple,  isolated  traps  such  as  T2  (cf  Fig.  3).  However,  it  rises  much  too  fast  to  explain  the  capture  process  of 
the  pore-type  traps  (To,  in  Fig.  3).  The  logarithmic  approximation,  on  the  other  hand,  works  fairly  well  for 
filling  fractions  up  to  about  0.5,  but  fails  beyond  that  point.  From  the  exact  solution,  it  is  seen  that  even  at  tp 
=  20  ms,  complete  saturation  has  not  taken  place.  Also,  in  emission,  the  exact  solution  is  much  slower  than 
the  exponential  solution,  because  at  higher  values  of  f  the  emitting  electrons  experience  a  strong  coulomb 
barrier  and  are  slowed  down.  In  Fig.  2,  we  have  also  simulated  a  boxcar  analysis  on  the  BioRad  DL4600 
instrument  by  indicating  a  common  set  of  sampling  points,  t|  =  61.0  ms  and  t2  -  152.6  ms,  referenced  with 
respect  to  tp.  This  choice  leads  to  an  emission  rate  of  In(t2/t|)/(t2  -  ti)  =  10  s’’  at  the  signal  maximum  of  a  trap 
such  as  T2  (Fig.  3),  which  has  an  exponential  emission  (Eqs.  8  and  9).  However,  the  emission  for  trap  To  is 
far  from  exponential,  so  that  the  “standard”  analysis  will  be  highly  inaccurate  in  this  case. 

The  experimental  DLTS  data,  for  filling  pulse  lengths  of  0.2, 1.0,  5.0,  and  20  ms,  are  shown  as  dashed 
lines  in  Fig.  3.  Here  we  have  plotted  AC/C,  where  AC  =  C(ti)  -  C(t2),  and  C  is  the  equilibrium  capacitance 
under  the  reverse-bias  condition,  Vr  =  -5  V.  It  can  be  shown  that  -AC/C  =  F;;^,Nt/2Nd,  where  Nj  is  the  trap 
concentration,  and  F^.  is  a  factor  which  is  close  to  unity  for  small  trap  concentrations  (Nt  «  No)  and 
energies  that  are  not  too  deep,  but  <  1  otherwise  [3].  For  our  case,  Nt  =  4Ttrp^NssNp,  where  Np  is  the  volume 
density  of  pores.  From  the  TEM  measurements,  the  sheet  density  of  pores  is  about  3  x  10'°  cm'^,  and  the 
volume  density  Np  is  then,  very  approximately,  (3x10)  si5xl0  cm'  .  Thus,  Nt  «  6  x  10  cm  ,  and 
from  this  value  and  also  Et  =  Ess  -  0.8  eV,  we  can  calculate  F^  =  0.25  [3].  The  actual  DLTS  signal  is  oc 
NT[f(ti)  -  f(t2)],  as  shown  in  Fig.  2,  and  f(ti)  -  f(t2)  is  calculated  to  be  0.646  at  the  peak  of  the  20-ms 
theoretical  curve,  in  Fig.  3.  Thus,  from  the  TEM  data,  we  would  predict  that  FxNT[f(ti)  -  f(t2)]/ND  «  0.10, 
whereas  we  actually  need  a  value  of  about  0.03  to  fit  the  data  at  the  peak,  as  shown.  In  other  words,  we  need 
an  Nt  value  of  about  2  x  lO’^  cm'^  to  fit  the  data,  which  is  not  outside  the  error  of  that  determined  by  TEM 
(6x10  cm' ),  considering  that  the  latter  value  is  a  rather  crude  estimate. 

The  normalization  factor  for  the  20-ms  curve  is  now  applied  to  the  other  three  theoretical  curves,  and 
they  reproduce  their  respective  experimental  peak  magnitudes  quite  well.  Furthermore,  the  temperature 
shifts  are  also  well  reproduced,  giving  strong  validity  to  our  model.  Finally,  both  the  experimental  and 
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theoretical  curves  become  narrower  at  larger  values  of  tp.  The  experimental  curves  are  of  course  broader 
than  their  theoretical  counterparts,  because  we  have  not  considered  the  known  variations  in  rp,  and  the 
possible  variations  in  Ess-  The  variation  in  rp  is  not  likely  the  cause  of  the  broadening,  because  neither  a 
doubling  nor  a  halving  of  the  pore  size  moves  the  curves  by  more  than  a  few  K.  On  the  other  hand,  an 
increase  of  Ess  from  0.80  to  0.85  eV  moves  the  curves  up  by  almost  20  K,  which  is  sufficient  to  explain  the 
line  broadening.  Indeed,  a  ±0.05  eV  variation  in  Ess  seems  quite  reasonable,  since  some  of  the  pores  will 
undoubtedly  be  close  enough  to  influence  each  other.  It  also  should  be  noted  that  adding  more  traps  at 
different  values  of  Ess  would  also  bring  the  total,  fitted  Nj  closer  to  the  TEM  estimate. 

2.4  Microcracks  in  GaN 

Freestanding  layers  of  GaN,  grown  on  AI2O3  by  hydride  vapor  phase  epitaxy  and  then  separated 
from  the  AI2O3,  are  beginning  to  be  used  as  substrates  for  GaN-based  devices.  Thus,  defects  in  these  layers 
are  of  both  scientific  and  practical  interest.  Because  these  layers  are  thick,  they  are  subject  to  cracking,  and 
this  problem  must  be  solved  before  they  can  become  useful.  Recently,  we  have  imaged  freestanding  GaN 
with  scanning  surface  potential  microscopy  (SSPM)  [12],  and  have  found  features  that  resemble  subsurface 
microcracks,  in  that  they  project  onto  the  surface  as  a  straight,  narrow  strip,  of  width  about  2  pm.  The  SSPM 
image  of  a  40-pm  x  40-pm  area  is  shown  in  Fig.  4,  and  a  line  scan  of  the  potential,  in  Fig.  5.  Note  that  this 
potential  line  scan  has  a  shape  much  like  the  energy  line  scan  illustrated  in  Fig.  1  (except  inverted).  We  can 
attempt  to  model  the  microcrack  as  a  thin,  infinite-area  box  with  charge  on  either  side,  of  sheet  density  Nss- 
Since  charge  balance  requires  that  Nss  =  Nqw,  and  since  w  «  2  pm,  from  the  figure,  and  No  «  n  «  6  x  lO'^ 
cm'^,  from  Hall-effect  measurements,  we  find  that  Nss  »  1  x  lO'^  cm'^.  Then,  the  potential  at  the  surface 
should  be  ({)  »  eNssV2sND  »  16  V.  This  value  is,  of  course,  much  higher  than  what  is  actually  measured, 
possibly  because  the  measurement  takes  place  at  the  “end”  of  the  box,  at  the  surface,  and  here  the  potential 
may  be  much  reduced  from  the  maximum  value,  say,  that  toward  the  middle  of  the  box.  Thus,  if  s  difficult 
to  be  quantitative  in  this  case,  but  nevertheless,  the  line  scan  shown  in  Fig.  5  clearly  is  giving  information  on 
the  lateral  shape  of  the  potential  around  a  crack.  To  be  fully  understood,  this  type  of  problem  will  have  to  be 
explored  in  more  detail. 

2.5  Dislocations  in  GaN 
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Chems  and  coworkers  have  developed  electron  holography  in  conjunction  with  a  transmission 
electron  microscope  to  study  charge  on  threading  dislocations  [13,14],  Basically,  the  interference  fiinges  in 
the  hologram  are  affected  by  an  electric  field  in  the  vicinity  of  the  sample,  and  this  field  is  of  course  directly 
related  to  the  potential.  The  potential  of  a  cylindrical  charge,  given  by  Eq.  3,  is  directly  applicable  to  an 
open-core  screw  dislocation,  of  radius  rp,  with  an  acceptor  sheet  density  of  Nss  on  the  inner  surface. 
However,  to  model  an  edge  dislocation  with  this  formula,  we  must  imagine  a  small  open  core  of  radius  rp, 
and  then  define  a  line  charge  as  e27rrpNssL,  where  L  is  the  dislocation  length.  From  end-on  images  of  typical 
threading-edge  dislocations  in  GaN,  it  seems  that  a  good  value  for  rp  is  about  2  A.  The  important  parameter 
for  comparison  with  theory  and  the  results  of  other  experiments  is  the  charge  per  unit  length,  e2TcrpNss.  Even 
more  eonveniently,  it  can  be  normalized  to  the  c-axis  parameter,  as  e27rrpNssc,  where  c  =  5.185  A.  Chems’ 
data  (Fig.  3  of  Ref  14)  are  presented  in  Fig.  6,  and  our  fit  of  these  data  is  also  shown.  For  this  sample.  No  ~ 
lO'^  cm'^,  and  the  best  fit  of  the  potential  is  then  found  for  Nss  ~  2.6  x  lO'"^  cm'^,  giving  a  normalized  line 
charge  e2TrrpNssc  =  1 .7  e  per  c-lattice  distance  (e/c).  Chems  also  attempted  a  fit,  but  with  a  much  simpler 
model  of  the  potential,  neglecting  screening.  Nevertheless,  his  fitted  line  charge  was  about  2e/c,  not  much 
different  than  ours.  It  should  also  be  mentioned  that  we  have  added  a  positive  background  potential  of  2.6  V 
in  order  to  fit  the  data  more  accurately.  It  would  be  interesting  to  apply  a  pulsed  external  potential  to  this 
dislocation  and  study  the  dynamics  of  carrier  trapping  and  emission,  as  was  done  for  the  pores  in  SiC, 
considered  earlier. 
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Fig.l  Conduction  band  profiles  under  reverse  bias,  and  forward  bias,  respectively.  The  dashed  lines  denote 
the  bands  after  the  filling  pulse  (curve  2,  at  the  end  of  the  filling  pulse  but  still  in  forward  bias;  and  curve  3, 
immediately  after  reapplication  of  the  reverse  bias). 


Fig.  2  Fractional  occupation,  during  capture  and  emission  at  350  K,  of  a  pore  in  porous  SiC.  An  exact 
calculation  is  compared  with  exponential  and  logarithmic  approximations.  The  filling  pulse  length  is  tp,  and 
the  sampling  points  on  the  emission  transient  are  ti  and  t2,  respectively. 
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Fig.  3  Experimental  (dashed  lines)  and  theoretical  (solid  lines)  DLTS  curves  for  different  filling-pulse 
lengths,  0.2,  2,  5,  and  20  ms,  in  porous  SiC.  Trap  T2  is  a  “normal”  trap  (impurity  or  point  defect),  which 
obeys  exponential  kinetics,  and  trap  To  is  related  to  the  pores. 


Fig.  4  AFM  and  SSPM  images  of  microcrack-like  defects  in  HVPE  GaN.  A  line  scan  of  the  surface 
potential,  shown  in  Fig.  5,  has  a  path  denoted  by  the  dotted  line  in  this  figure. 
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Fig.  5  Lateral  profile  of  the  surface  potential  along  the  dotted  line  shown  in  Fig.  4.  Note  that  the  line  crosses 
two  cracks. 


Distance  from  dislocation  core  (nm) 

Fig.  6  Lateral  potential  profile  across  a  threading-edge  dislocation  in  GaN.  The  dotted  line  is  a  theoretical 
fit,  assuming  a  line  charge  of  1.7  e/c. 
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3.0  Characterization  of  bulk-like  GaN  grown  by  hydride  vapor  phase  epitaxy 

3.1  Introduction 

GaN  grown  by  hydride  vapor  phase  epitaxy  (HVPE)  has  a  high  growth  rate  (~  100  pm  per  hr),  and 
thus  is  capable  of  growing  thick  material,  more  than  1000  pm  in  some  cases  [1-5],  Although  these  layers  are 
usually  grown  on  AI2O3  substrates,  they  can  be  easily  separated  from  the  AI2O3  by  a  laser  irradiation 
technique  [4],  Recently,  separated  (or  unseparated)  HVPE  wafers  have  been  proposed  as  a  solution  to  the 
GaN  substrate  problem  because,  unlike  the  cases  in  Si,  GaAs,  SiC,  and  ZnO,  large-area  wafers  of  GaN 
cannot  be  obtained  by  the  usual  bulk-growth  techniques.  Thus,  it  is  important  to  characterize  the  electrical 
properties  of  these  HVPE  GaN  templates,  and,  in  particular,  to  identify  and  determine  the  concentrations  of 
their  dominant  donors,  acceptors,  and  traps.  It  is  well  known  that  HVPE  GaN  layers  grown  on  AI2O3  can 
often  be  described  in  terms  of  two  regions,  a  highly  conductive  layer  adjacent  to  the  GaN/Al203  interface,  of 
thickness  0.2  -  50  pm,  and  the  rest  of  the  sample,  comprising  the  “bulk”  portion,  of  much  lower 
conductivity  [6].  The  interface  layer  has  been  extensively  studied  in  the  past,  and  contains  very  high 
concentrations  of  On  donors  and  Ga- vacancy  acceptors,  as  shown  in  Fig.  1,  and  also  threading  dislocations 
[7-10].  However,  thick  layers  usually  have  a  strong  bow  after  growth,  and  have  to  be  ground  flat  before 
being  used  as  substrates  for  subsequent  growth  of  an  active  layer  [3,4].  This  grinding/polishing  process 
typically  removes  all  or  part  of  the  low-quality  interface  layer  so  that  only  the  higher-quality  bulk  region 
remains.  We  would  normally  like  the  bulk  region  to  be  as  pure  as  possible  in  order  to  avoid  diffusion  of 
impurities  and  defects  into  the  active  layer  and  the  trapping  of  carriers  injected  from  the  active  layer  into  the 
substrate.  It  is  the  purpose  of  this  work  to  review  our  present  knowledge  of  donors,  acceptors,  and  traps  in 
thick  HVPE  GaN  and  present  some  new  results  and  interpretations.  We  will  deal  principally  with  HVPE 
GaN  layers  grown  at  three  different  laboratories,  designated  as  A,  B,  and  C.  Lab-A  and  Lab-C  samples  are 
1-pm  and  30-pm  thick,  respectively,  and  are  still  on  their  AI2O3  substrates.  The  Lab-B  sample,  on  the  other 
hand,  was  originally  500-pm  thick,  and  then  was  separated  from  its  AI2O3  substrate  and  ground  and 
polished  flat  to  a  thickness  of 250  pm. 

3.2  Secondary  ion  mass  spectroscopy:  donor-type  impurities 

Secondary-ion  mass  spectroscopy  (SIMS)  measurements  [10]  of  [Si],  [O],  and  [C]  have  been  carried 
out  on  the  Lab-A  and  Lab-B  samples.  In  the  bulk  regions,  both  samples  show  [O]  and  [C]  in  the  mid  lO'^ 
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cm'^  range,  but  the  Lab-A  sample  has  [Si]  >  [O]  in  this  region  (see  Fig.  1),  and  the  Lab-B  sample,  [O]  >  [Si]. 
Another  group  has  carried  out  SIMS  measurements  on  material  from  Lab  B,  and  has  also  found  that  [O]  > 
[Si].  From  temperature-dependent  Hall-effect  (T-Hall)  measurements,  discussed  below,  we  find  that  the 
total  shallow  donor  concentration  in  the  Lab-B  sample  is  only  about  8x10  cm' ,  a  lower  value  than  [0]  + 
[Si],  as  determined  by  SIMS.  However,  SIMS  measurements  are  not  always  accurate  at  these  low 
concentrations,  and,  also,  some  of  the  O  and  Si  may  not  be  electrically  active  as  donors.  The  other  common 
impurity,  C,  would  be  expected  to  substitute  for  N  and  be  an  acceptor  in  n-type  GaN,  but  the  SIMS  value  of 
[C],  mid-lO'^  cm'^,  is  much  greater  than  the  acceptor  concentration  Na,  low-lO'^  cm'^,  determined  by  the 
Hall  effect;  thus,  either  the  SIMS  data  are  not  accurate  due  to  background  limitations,  or  most  of  the  C  is  not 
electrically  active.  In  this  regard,  it  should  be  pointed  out  that  the  profiles  of  C  and  O  in  some  of  these 
samples  are  quite  similar,  suggesting  a  correlation  between  the  two,  either  in  the  sample  or  the  background. 

3.3  Positron  annihilation  spectroscopy  (PAS):  acceptor-type  defects 

Positrons  injected  into  defect-free  GaN  are  annihilated  by  electrons  in  a  mean  time  of  160  -  165  ps. 
However,  if  there  are  negatively  charged  vacancies  present,  some  of  the  positrons  will  become  trapped  at 
these  locations,  and  will  have  longer  lifetimes,  because  of  the  reduced  electron  density  at  vacancies.  In  the 
case  of  GaN,  the  Ga  vacancy  Voa  (but  not  the  N  vacancy)  would  be  expected  to  fill  this  role,  and  indeed, 
PAS  has  been  used  to  identify  and  quantify  Voa-related  defects  [9].  In  fact,  comparisons  of  Voa 
concentrations  with  acceptor  concentrations  Na  in  a  series  of  Lab-A  HVPE  GaN  samples,  with  Na  ranging 
from  lO’^  to  lO’^  cm'^,  show  that  [Voa]  »  Na,  to  within  experimental  error  [7,9].  Also,  a  sample  from  Lab  B 
has  been  shown  to  have  [Vca]  «  2  x  10*^  cm'^  [9],  very  close  to  our  value  of  Na  determined  by  T-Hall 
measurements,  discussed  below.  Thus,  it  appears  that  Voa,  and  not  any  impurity,  is  the  dominant  acceptor  in 
HVPE  GaN.  Indeed,  theory  predicts  that  Voa  centers  should  be  abundant  in  n-type  GaN  [11].  However,  in 
GaN  grown  by  some  other  techniques,  such  as  by  low-pressure  metal-organic  vapor-phase  epitaxy,  the 
dominant  acceptor  is  Cn,  rather  than  Vca- 

3.4  Photoluminescence  -  donors 

A  4-K  PL  spectrum  of  the  near-band-edge  (exciton)  region,  3.43  -  3.48  eV,  is  shown  in  Fig.  2  for  a 
HVPE  GaN  sample  from  Lab  B.  The  sharp  lines  at  3.47123,  3.47225,  and  3.47305  eV  are  likely  neutral 
donor-bound  A  excitons  (D'^Xa’s),  while  the  broader  line  at  3.4792  eV  is  the  free  A  exciton  Xa.  The  line  at 
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3.47609  eV  may  be  an  excited  (rotator)  state  of  a  D®Xa  species,  or  possibly  a  D^Xb  transition.  It  has  been 
reported  that  the  D°Xa  line  in  unstrained  material  should  lie  at  3.471  -  3.472  eV,  a  result  that  suggests 
immediately  that  the  present  Lab-B  sample  does  not  have  a  high  strain.  This  is  expected,  since  strain 
decreases  with  thickness,  and  this  sample  is  250-pm  thick.  The  full  width  at  half  maximum  (FWHM)  for 
each  of  these  D°X  lines  is  about  0.4  meV,  which  indicates  excellent  quality.  Another  group  of  PL  lines 
appears  in  the  region  3.44  -  3.46  eV,  with  a  strong  line  at  3.44686  eV  and  a  somewhat  weaker  one  at  a 
slightly  higher  energy,  3.44792  eV.  Most  workers  agree  that  the  3.44686-eV  line  is  a  two-electron  satellite 
(TES)  replica  of  one  of  the  three  D°X  transitions,  although  there  is  disagreement  over  which  one.  The 
3.44792-eV  line  would  then  be  a  TES  transition  from  an  excited  (rotator)  state  of  that  particular  D^’X  line. 
However,  there  are  also  a  few  weaker  lines  at  even  lower  energies,  and  two  of  these  are  designated  n=3  and 
n=4  in  Fig.  2.  We  have  earlier  pointed  out  [12]  that  if  we  combine  these  two  lines  with  the  strongest  n=l 
D°X  line  at  3.47225  eV,  and  the  strongest  n=2  line  at  3.44686,  then,  as  seen  in  the  inset  of  Fig.  2,  all  four 
peak  energies  closely  follow  the  hydrogenic  relationship,  EpL,n  =  Epl,i  -  Rexpt(l  -  1/n^)  =  3.47225  - 
0.0339(l-l/n^),  so  that  Rexpt  =  33.9  meV.  Rexpt  is  the  “experimental”  Rydberg  for  GaN.  The  theoretical 
Rydberg  is  derived  from  effective-mass  theory:  Ruico  =  13.6m*/so^  eV,  where  m*  is  the  electron  effective 
mass  relative  to  the  rest  mass,  and  so  is  the  dielectric  constant  relative  to  the  value  in  vacuum.  The  presently 
accepted  values  of  m*  and  so  are  0.22  [13]  and  9.4  [14],  respectively,  giving  Rtheo  =  33.9  meV,  fortuitously 
close  to  Rexpt-  The  clear  and  simple  interpretation  of  these  results  is  that  the  donor  represented  by  the  D'^X 
line  at  3.47225  eV  is  truly  hydrogenic,  i.e.,  effective-mass-like.  From  the  SIMS  results,  the  two  obvious 
candidates  for  this  donor  are  Sica  and  On;  however,  it  would  seem  that  On  is  by  far  the  better  of  the  two 
candidates,  because  O  and  N  have  nearly  equal  ionic  radii,  and  thus  only  a  small  lattice  distortion  would  be 
expected  from  the  replacement  of  N  with  O.  In  other  words,  the  potential  in  the  vicinity  of  On  should  be  just 
that  of  the  original  Nn  plus  a  point  positive  charge  at  the  core,  which  should  lead  to  a  hydrogen-like 
potential.  Interestingly,  a  recent  calculation  based  on  effective-mass  theory  [15]  shows  that  the  details  of  the 
band-edge  wave  function  don’t  make  much  difference  for  the  On  donor  energy,  whereas  they  do  for  the 
energies  of  other  donors,  such  as  Sioa-  This  fact  would  seem  to  support  the  idea  that  On  is  hydrogenic  or 
nearly  so.  In  this  same  paper,  the  hydrogenic  ground  state  energy  (Rydberg)  was  given  as  Rtheo  =  29.6  meV; 
that  number  was  evidently  derived  from  the  usual  formula  (13.6mVso^  ),  but  with  older  values  of  dielectric 
constant  so  =  9.8,  and  effective  masses  nv=  0.19  and  mx  =  0.22  [giving  an  average  m*  =  (nvmx^)’^^  = 
0.209].  Our  values  of  m*  =  0.22  and  so  =  9.4  both  come  from  recent  experiments,  and  thus  our  value  Ruieo  = 
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33.9  meV  may  be  more  accurate.  In  any  case,  the  origin  of  the  difference  between  the  values  29.6  and  33.9 
meV  is  well  understood.  Another  disagreement,  however,  which  is  not  understood,  involves  a  comparison 
of  2p  and  3p  IR-absorption  states  [13],  which  leads  to  an  experimental  value  of  the  Rydberg,  Rexpt  -  29.1  ± 
0.5  meV,  much  smaller  than  our  value,  Rexpt  =  33.9  meV.  This  difference  will  require  further  investigation. 

Another  controversial  issue  is  concerned  with  the  identification  of  the  dominant  DOX  line  at  3.47225 
eV.  It  was  originally  identified  as  On  rather  than  Sica  mainly  on  the  basis  of  SIMS  data,  i.e.,  because  these 
data  show  more  O  than  Si  in  nearly  all  samples  from  Lab  B  [16].  However,  as  mentioned  earlier,  a  SIMS 
measurement  of  a  particular  element,  say  O,  includes  all  of  that  element,  not  just  the  substitutional  fraction 
(On  in  this  case).  In  this  regard,  we  note  that  the  SIMS  values  of  [O]  and  [Si]  in  Lab-B  samples  are  always 
significantly  higher  than  the  total  donor  concentration,  as  measured  by  T-Hall  measurements.  Thus,  SIMS 
measurements  are  not  conclusive  here.  To  further  complicate  the  O-vs.-Si  matter,  another  group  [17]  has 
compared  the  transient  behaviors  of  the  3.47225-eV  (n=l)  and  3.44686-eV  (n=2)  lines,  and  has  concluded 
that  these  lines  cannot  have  a  common  origin.  Instead,  they  associate  the  3.47 1 23 -eV  (n=l)  line  with  the 
3.44686-eV  (n=2)  line,  and  assign  both  to  On.  But  a  brief  glance  at  Fig.  2  illustrates  a  potential  problem 
with  this  assignment,  namely  that  the  3.47123-eV  (n=l)  line  is  by  far  the  weakest  D°X  line,  and  the 
3.44686-eV  (n=2)  line,  by  far  the  strongest  TES  line.  In  the  present  work,  we  are  suggesting  that  the 
3.47225-eV  (n=l),  3.44686-eV  (n=2),  3.44215-eV  (n=3),  and  3.44051-eV  (n=4)  lines  all  belong  to  On,  not 
on  the  basis  of  SIMS  or  transient  PL,  but  because  these  four  lines  seem  to  fit  a  33.9-meV  hydrogenic  model 
very  well,  and  that  the  best  candidate  for  a  hydrogenic  donor  is  On,  rather  than  Sioa-  However,  at  this  point, 
we  cannot  explain  the  experimental  Rydberg  of  29.1  meV,  derived  from  the  2p  and  3p  absorption  lines  [13], 
or  the  difference  in  transient  PL  behavior  of  the  3.47225-eV  (n=l)  and  3.44686-eV  (n=2)  PL  lines  [17]. 
Clearly,  more  work  needs  to  be  done  to  clear  up  these  matters. 

3.5  Hall-effect  measurements:  donors  and  acceptors 

In  a  simple  picture,  conductivity  is  defined  by  a  =  nepc,  where  n  is  the  carrier  (electron) 
concentration  (not  to  be  confused  with  the  symbol  for  quantum  number,  used  previously),  and  pc  is  the 
conductivity  mobility,  defined  by  pc  =  e(T)/m*.  Here,  x  is  a  “relaxation  time”,  obeying  [18] 

x\E)  =  X^;\E)  +  Xpe''(E)  +  Xpo''(F: )  +  +  Tdis‘‘(£'  )  (1) 
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and  (x)  denotes  an  average  over  energy  E.  The  relaxation  time,  x(E),  depends  upon  how  the  electrons 
interact  with  the  lattice  vibrations  as  well  as  with  extrinsic  elements,  such  as  charged  impurities  and  defects. 
For  example,  acoustical-mode  lattice  vibrations  scatter  electrons  through  the  deformation  potential  (xac)  and 
piezoelectric  potential  (Xpe);  optical-mode  vibrations  through  the  polar  potential  (Xpo);  ionized  impurities  and 
defects  through  the  screened  coulomb  potential  (xjj);  and  charged  dislocations,  also  through  the  coulomb 
potential  (xdis).  The  strengths  of  these  various  scattering  mechanisms  depend  upon  certain  lattice 
parameters,  such  as  dielectric  constants  and  deformation  potentials,  and  extrinsic  factors,  such  as  donor, 
acceptor,  and  dislocation  concentrations.  No,  Na,  and  Njis,  respectively  [18-20].  In  reality,  a  closed-form 
expression  for  Xpo  cannot  be  written,  but  suitable  approximations  are  available  [21]. 

In  general,  it  is  difficult  to  measure  pc,  but  much  easier  to  measure  pn  =  Rh<7,  where  Rh  is  the  Hall 
coefficient.  In  the  relaxation-time  approximation,  ph  is  calculated  from  pH  =  e(x^)/m*(x)  =  rpc.  It  is  also 
convenient  to  define  a  “Hall  concentration”,  nn  =  n/r  =  -  I/cRh.  Thus,  a  combined  Hall-effect  and 
conductivity  measurement  gives  nn  and  ph,  although  we  would  prefer  to  know  n,  not  nn;  fortunately, 
however,  r  is  usually  within  20%  of  unity,  and  is  almost  never  as  large  as  two.  In  any  case,  r  can  often  be 
calculated  or  measured  so  that  an  accurate  value  of  n  can  usually  be  determined. 

The  fitting  of  pn  vs.  T  data  should  be  carried  out  in  conjunction  with  the  fitting  of  n  vs.  T  data,  and  the 
relevant  expression  here  is  the  charge-balance  equation  [18]: 


n  +  N^  = 


No 
l  +  n/ 


(2) 


where  we  have  assumed  only  one  type  of  donor,  with  a  single  charge  state,  and  where 


^DO 

(t>o=—e'‘  NcT^'^i 

S\ 


(3) 


Here,  go/gi  is  a  degeneracy  factor,  Nc'  =  2(27rmn*k)^^W  is  the  effective  conduction-band  density  of  states  at 
IK,  h  is  Planck’s  constant,  Ed  is  the  donor  ground-state  energy,  and  Edo  and  ttD  are  defined  by  Ed  =  Edo  - 
ttDT.  If  more  than  one  donor  is  needed  to  fit  the  data,  then  equivalent  terms  are  added  on  the  right  hand  side 
of  Eq.  2.  As  mentioned  earlier,  examples  of  common,  single-charge-state  donors  in  GaN  are  Sioa  and  On.  If 
there  are  double  or  triple  donors,  or  more  than  one  acceptor,  proper  variations  of  Eq.  2  can  be  found  in  the 
literature  [18]. 
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If  the  donors  are  effective-mass-like,  they  will  have  a  set  of  excited  states,  much  like  those  of  hydrogen. 
Using  standard  statistical  theory,  we  can  add  hydrogenic-type  excited  states  (j  =  2,  3,....,m)  to  the  analysis 
by  modifying  (j)D  [18,22]. 


^Z) 


=  ^e  * 
^1 


^DO 


(4) 


where  we  have  assumed  that  gj/go  =  f,  as  is  the  case  for  the  hydrogen  atom,  and  also  that  ap  is  the  same  for 
each  state.  (Actually,  in  any  case,  ao  should  be  small  for  an  effective-mass-like  donor  state.)  At  low 
temperatures,  only  the  ground  state  will  be  occupied,  and  the  additional  term  in  the  denominator  of  Eq.  4 
will  be  small.  However,  at  higher  temperatures,  the  n  vs.  T  curve  will  be  modified  [12]. 

The  Lab-B  Hall  mobility  pn  data,  and  the  theoretical  fit  are  plotted  vs.  temperature  in  Fig.  3.  From  these 
data,  an  acceptor  concentration  Na  =  1.3  x  lO'^  cm‘^  is  deduced,  the  lowest  ever  determined  in  GaN.  The 
carrier  concentration  data,  corrected  for  the  Hall  r-factor,  are  plotted  in  Fig.  4,  along  with  the  theoretical  fit 
(Eq.  2).  Here  two  donors  are  found  from  the  fit:  Ndi  =  7.8  x  lO'^  cm’^,  Edi  =  28.1  meV;  andND2  =  1.1  x  lO’^ 
cm‘^,  Ed2  =  53.2  meV.  Also,  the  fitted  acceptor  concentration  Na  is  7.2  x  lO’"*  cm'^,  a  little  smaller  than  the 
value  found  from  the  mobility  fit,  but  not  considered  to  be  as  accurate  as  the  latter.  Note  that  ED(Hall)  < 
Ed(PL)  by  a  few  meV,  as  is  expected  from  the  effects  of  waveflinction  overlap  and  excited  donor  states,  as 
discussed  in  detail  elsewhere  [12]. 

The  most  straightforward  interpretation  of  all  of  the  above  results,  is  that  the  Lab-B  sample  contains  On 
at  a  concentration  of  8  x  lO'^  cm'^,  and  that  On  forms  a  hydrogenic  donor  of  energy  33.9  meV  from  the 
conduetion  band,  with  distinctive  PL  lines  and  a  Hall-effect  energy  a  little  smaller,  28.1  meV,  in  this 
particular  case. 
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Depth  (|im) 

Fig.  1.  The  interface  region  of  a  1.3-|j.m-thick  Lab-A  HVPE  GaN  layer  grown  on  an  AI2O3  substrate.  Shown 
are  SIMS  profiles  for  O  and  Si,  an  electrochemical  C-V  profile  for  the  electron  concentration  n,  and  the 
average  Ga-vacancy  Voa  concentration  deduced  from  positron  annihilation  spectroscopy. 


Fig.  2.  Photoluminescence  spectrum  in  the  donor-bound-exciton  (D®X)  and  two-electron  satellite  regions. 
The  four  lines  marked  D°X  n=l,  etc,  have  been  plotted  in  the  inset  and  obey  the  relationship  EpL,n  =  3.47225 
-0.0339(l-l/n^). 
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Fig.  4.  Carrier  concentration  vs.  inverse  temperature  for  a  Lab-B  HVPE  GaN  layer.  The  solid  line  is  a 
theoretical  fit. 
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Fig.  5.  A  comparison  of  DLTS  data  for  two  HVPE  samples,  and  an  RMBE  sample. 
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Fig.  6.  Effects  of  electron  irradiation  (1  MeV,  1  x  lO'^  cm'^)  on  HVPE-C  GaN/Al203  layer. 
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